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The nuclei 163164 Dy have been investigated using the Oslo method on data from the pick-up reaction 
164 Dy( 3 He, ay) 163 Dy and the inelastic scattering I64 Dy( 3 He, 3 He'y) 164 Dy, respectively. The radiative strength 
functions for both nuclei have been extracted, and a small resonance centered around Ey ~ 3 MeV is observed in 
both cases. The parameters of this so-called pygmy Ml resonance (the scissors mode) are compared to previous 
results on 160 > 161 ' 162 Dy using the Oslo method, and to data on 163 Dy measured by the Prague group using the 
two-step cascade method. In particular, the integrated reduced transition probability 6(M1 f) of the pygmy res- 
onance is compared with neighboring dysprosium isotopes. We also observe an enhanced strength in the region 
above Ey ~ 5 MeV in 164 Dy. Possible origins of this feature are discussed. 

PACS numbers: 25.20.Lj, 24.30.Gd, 25.55.Hp, 27.70.+q 



I. INTRODUCTION 

A continuing effort has long been devoted to study y decay 
from excited nuclei. The radiative strength function (RSF) 
represents the mean value of the decay probability via a 7 ray 
with energy Ey, and contains rich information on the average 
electromagnetic properties of the nucleus. 

For high-energy 7 transitions (^7 — 20 MeV), the RSF is 
dominated by the giant electric dipole resonance (GEDR). 
At lower energies other resonances have been discovered, 
such as the giant magnetic dipole resonance (GMDR, also 
called the giant magnetic spin-flip resonance) and the elec- 
tric quadrupole resonance; however, these have a significantly 
lower strength (T). In addition, there are other structures ob- 
served in the RSF governed by various collective modes of 
the nucleus. These structures are often referred to as pygmy 
resonances due to their low strength compared to the GEDR. 
There are two known pygmy resonances: the El resonance 
for 7 ray energies between 5—10 MeV that is believed to 
stem from neutron skin oscillation [2|, and the Ml resonance 
called the scissors mode, which is observed in the region of 
Ey = 3 MeV for rare-earth nuclei 0. 

The RSFs below the neutron threshold have been studied 
mainly by (7, /) experiments, also called nuclear resonance 
fluorescence (NRF) [4|. Other methods, such as the two-step- 
cascade (TSC) method [5 1 and the Oslo method (6), have also 
successfully provided data on the RSF for many nuclei. The 
latter method enables us to extract the RSF for 7 ray energies 
up to the neutron binding energy B„. This method has been 
used in the present analysis. 

Previously, experiments have been performed on 
160,161, i62 Dy J7J using the Oslo method ©. From these 
data, the widths of the Ml pygmy resonance have been found 
to be about two times greater than the width found for 163 Dy 
obtained by the Prague group using the TSC method JS). In 
the present work we have studied 163 164 Dy to investigate the 
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discrepancy between the measured widths. In particular, we 
have compared the total integrated strength B(M1 f) of all the 
mentioned Dy isotopes. 

Details about the experimental method are presented in 
Section II, followed by the experimental results for the RSF 
in Section III. Finally, conclusions will be drawn in Section 
IV. 



II. EXPERIMENTAL PROCEDURE AND DATA ANALYSIS 

The experiment was conducted at the Oslo Cyclotron Lab- 
oratory (OCL), using a 38-MeV beam of 3 He particles. 
The target of 98.5% enriched 164 Dy had a thickness of 
1.73 mg/cm 2 , and the reactions 164 Dy( 3 He, G!7) 163 Dy and 
164 Dy( 3 He, 3 He'7) 164 Dy were studied. 

The 7 rays and ejectiles were measured with the CACTUS 
multidetector array (5), which consists of a sphere of 28 col- 
limated Nal 7 detectors with total efficiency of 15% of An, 
surrounding a vacuum chamber which contains eight AE — E 
Si particle telescopes with thicknesses of 140 and 1500 jim, 
respectively. The particle telescopes were placed in forward 
direction at 45° relative to the beam axis. 

From the known Q-values the excitation energy of the nu- 
clei were calculated from the detected ejectile energy using 
reaction kinematics. The particles and 7 rays were measured 
in coincidence; hence, each 7 ray can be assigned to an ini- 
tial excitation energy of the nucleus. The 7 ray spectra were 
unfolded using the known response functions of the CACTUS 
detector array iflOl . The excited nuclei will decay through a 
cascade of 7 rays down to the ground state. By using the first- 
generation method ATI , we were able to isolate the first (pri- 
mary) 7rays emitted in each 7decay cascade. The distribution 
of primary 7 rays was found for each excitation energy bin, 
giving an excitation energy vs. 7 ray energy matrix denoted 
by P(Ei,Ey). The primary 7 ray spectrum was normalized to 
unity for each excitation energy bin, which then represents 
the decay probability for each 7 ray with energy Ey decaying 
from a certain excitation energy E,: J^' _ Emjn P(Ei,Ey) = 1. 

The primary 7 ray matrix for 164 Dy is shown in Fig. uj The 
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TABLE I: Parameters used for normalizing p and ST. 
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FIG. 1: The primary y ray matrix for 164 Dy, displaying primary y 
rays emitted at each initial excitation energy. 



diagonal line of the matrix corresponds to decay directly to 
the ground state (Ey = Ei); however, it is more probable to de- 
cay to the first excited 2 + and 4 + states. Therefore, a ridge is 
formed in the matrix for decay to £, ~ 150 keV, see Fig.[T] 

The original Brink- Axel hypothesis |[T2l ITJl l states that the 
GEDR can be built on every excited state, and that the prop- 
erties of the GEDR do not depend on the temperature of the 
nuclear state on which it is built. This hypothesis can be gen- 
eralized to include any type of collective excitation. Provided 
that this hypothesis is valid, the primary y ray matrix can be 
factorized as 



P(E i ,E r )cc^(E y )p(E i -E Y ), 



(1) 



where P(Ei,Ey) is the experimentally obtained and normal- 
ized primary y ray matrix. The function 2?(Ey) represents the 
radiative transmission coefficient and p (Ei — Ey) is the level 
density at the final energy Ef = Ej — Ey. The above factor- 
ization is based on the essential assumption that the system 
is fully thermalized prior to y emission, so that the reaction 
can be described as a two-stage process of which the first is 
the formation of the compound nucleus, which subsequently 
decays in a manner that is independent of the mode of forma- 
tion 04). The formation of a complete compound state is as 
fast as ~ 10~ 18 s, significantly less than the typical life time 
of a state in the quasi-continuum which is ~ 10~ 15 s. There- 
fore, the assumption is believed to be reasonable, and the de- 
cay process is at least mainly statistical. Recently, it has been 
shown that Eq. ([T]l can be valid even in some cases where full 
thermalization is not achieved |[T5ll . 



However, there is experimental evidence that the Brink- 
Axel hypothesis is violated for high temperatures (above 1-2 
MeV). In particular, the width of the GEDR has been shown 
to depend on the temperature of the final states 11161 . For our 
experimental conditions the excitation energy (and thus the 
temperature) is relatively low and changes slowly with excita- 
tion energy (T ~ y / £/). Therefore, we assume that the radia- 
tive strength function does not depend on temperature in the 
energy region under consideration. 

The right hand side of Eq. ([T]) is normalized to unity, yield- 
ing 



P(Ei,E Y ) 



&{Ey)p(Ei-Ey 



..&{E')p(Ei-E' 



(2) 



Using this equation and applying a least squares fit to the pri- 
mary y ray matrix, a unique functional form of p (£, — Ey) 
and 3~(Ey) is derived [6 |, while the normalization is yet to be 
determined. There are infinitely many normalization options 
that reproduce the experimental primary y ray matrix. All the 
solutions are related to each other through the two transfor- 
mations |6| 



and 



p(Ei — Ey) =A exp[a(£, - Ey)] p (E, - E y ) (3) 



2T(Ey) = B exp(aEy) &(E r ), (4) 



where A, B and a are constants representing the absolute value 
of p(Ei~Ey) and ,^(Ey), and the slope of the two functions, 
respectively. These parameters are determined by normalizing 
Eqs. (|3]l and |4]) to known experimental data. The parameters 
A and a are identified by normalizing the experimental level 
density to known levels found from discrete spectroscopy at 
low energies. At higher excitation energies the experimen- 
tal level density is normalized to the level density determined 
from the known neutron resonance spacing data [ 1| at the neu- 
tron binding energy B„. The present experimental data extend 
up to about B„ — 1 MeV; an interpolation is thus required to 
reach B„. The back-shifted Fermi gas model [ 17 18 1 was ap- 
plied for this purpose: 



Pbs{E) = 7? 



exp(2y/aU) 



(5) 



where the constant 77 is applied to adjust Pb s (E) to the semi- 
experimental level density at B n . The intrinsic excitation en- 
ergy is given by U = E — C\— £p a i r , where C\ is the back-shift 
parameter equal to C\ = — 6.6A~° 32 MeV, where A represents 
the mass number. The pairing energy £p a ; r is based on the 
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FIG. 2: Upper panel: The level density of 164 Dy, normalized to 
known discrete levels and to p(B„) calculated from neutron res- 
onance spacing data, with an interpolation using the back-shifted 
Fermi gas model. Lower panel: the radiative transmission coeffi- 
cient of 164 Dy with extrapolations. The normalization is performed 
in the regions between the arrows. 



pairing gap parameters A p and A„ evaluated from odd-even 
mass differences [19] according to |20|. The parameter a = 
0.2 L4 - 87 MeV -1 corresponds to the level density parameter. 
The spin-cutoff parameter a is given by a 2 = 0.0888a7A 2 / 3 , 
where the nuclear temperature is described by 



y/uTc 



(6) 



The normalization of p (£, — Ey) for 1M Dy is displayed in the 
upper panel of Fig. [2] 

Finally, £F{Ey) is normalized by determining the coeffi- 
cient B, which gives the magnitude of 3?(Ey). We have the 
following relation between the total radiative width of neu- 
tron resonances (Fy) at the neutron binding energy and the 
radiative transmission coefficient ^(Ey) lETl : 



(r r ) = 



i 



\np{B„,J? 



J f 



£/ dEyB,y{Ey)p(B„ —Ey,Jf), 



(7) 

where D, ■ = 1 / p (B„ , /?) is the average spacing of .s-wave neu- 
tron resonances. The summation and integration extends over 
all final levels with spin Jf which are accessible by 7 radiation 
with energy Ey. 

Due to methodological difficulties, 3~{Ey) cannot 
be determined experimentally for low-energy 7 rays, 



statistics for 7 ray energies Ey > B„ — 1 MeV. We therefore 
extrapolate 2?(Ey) with an exponential function, as demon- 
strated for 164 Dy in the lower panel of Fig. [2] For further 
details of the normalization procedure, see Ref. [22]. The 
parameters used for normalizing p(Ei — Ey) and 9*(Ey) are 
given in Table [I] 

Note that the uncertainties displayed in Fig. [2] only reflect 
statistical uncertainties, and do not include the uncertainties 
related to the model used for normalization. This is also the 
case for the other figures showing experimental data. 



III. RADIATIVE STRENGTH FUNCTIONS 

Assuming that 7 decay taking place in the quasi-continuum 
is dominated by dipole transitions (L = 1), the radiative 
strength function can be calculated from the normalized trans- 
mission coefficient by 



f(E? 



1 ,J{Ey) 



In 



(8) 



Using the above relation, we obtain the experimental RSFs 
displayed in Fig. [3] We observe that they are increasing func- 
tions of 7 energy, and we can easily identify the Ml pygmy 
resonance in both cases. We expect the RFS to be composed 
of the pygmy resonance, the giant electric dipole resonance 
(GEDR) and the giant magnetic dipole resonance (GMDR). 
The Kadmenskii-Markushev-Furman (KMF) model [23] is 
employed to characterize the £1 strength. In this model, an 
excitation-energy dependence is introduced through the tem- 
perature of the final states 7V : 



fKMF 
J El 



(Ey) = 



1 



0.7 OEiEyTl^E* 



-An 2 T}) 



37t 2 h 2 c 2 



E E1 (E 2 -E 2 El ) 



(9) 



where CTgi, Tei and Eei denote the peak cross section, width 
and the centroid of the GEDR, respectively. In general, the 
KMF model describes experimental data very well; however, 
the temperature dependence violates the Brink- Axel hypothe- 
sis. In line with the previously mentioned argument that the 
temperature varies relatively little in our region of interest, we 
have assumed that the temperature can be considered to be 
constant. Thus, the Brink- Axel hypothesis is revived. 

It was found that a constant temperature of Tf — 0.3 MeV 
gives a good fit to the experimental data, in agreement with 
Ref. 0. For deformed nuclei the GEDR is split into two, and 
is therefore described as a sum of two strength functions given 
by Eq. ([9]). The GMDR is thought to be governed by the spin- 
flip Ml resonance 11221 . and can be described by a Lorentzian 
function: 



fMl{Ey) — 



1 



(TmiEvT 2 



r 1 mi 



3n 2 h 2 c 2 



E 2 " 



-E 2 Y 2 



(10) 



Ey < 1 MeV 1221 . In addition, the data suffer from poor 



where Om\, ^mi and Emi give the peak cross section, width 
and the centroid of the GMDR, respectively. The GEDR and 
GMDR parameters are taken from the systematics of Ref. HI 
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TABLE II: Parameters used for the radiative strength functions. 

Nucleus E El oj; r£[ £|j of, ffj £^ o^i Tmi ]hT 

(MeV) (mb) (MeV) (MeV) (mb) (MeV) (MeV) (mb) (MeV) 
lw Dy 12.37 278.50 3.17 15.90 139.04 5.12 iJi 1.49 4.00 0.300 
164 Dy 12.26 280.41 3.12 15.95 140.00 5.15 7.49 1.49 4.00 0.314 



TABLE III: Fitted pygmy-resonance parameters and normalization constants. 

NUCleUS Epy Gpy Tpy K 

(MeV) (mb) (MeV) 

lw Dy 2.81(9) 0.72(12) 0.86(19) 1.78(14) 
164 Dy 2.81(6) 0.53(6) 0.80(12) 1.72(6) 
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FIG. 3: Normalized RSFs of 163 164 Dy. 



calculated with the deformation parameter fa QJ. The Ml 
pygmy resonance f py is described by a Lorentzian function 
similar to the one given in Eq. ( |10) . All parameters are listed 
in Table II. 

The theoretical strength function is then given by 



f = K{f E i+fm)+f P . 



where fsi, /mi and f py represent the contributions from the 
GEDR, GMDR and the Ml pygmy resonance, respectively. 
The parameter K is a normalization constant. Together with 



the pygmy resonance parameters <j py , T py and 



E py , K is 



used as a free parameter when performing a least squares fit 
to adjust the total theoretical strength to the experimental data. 

The fit to the experimental data points is shown in Fig. [4] 
for both nuclei. The upper panels show the contributions 
KfEi and k/m\ and the sum of these two contributions. In 
the lower panels the sum k(/ei + /mi) is subtracted from the 
experimental data, and the fit to the Ml pygmy resonance 
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FIG. 4: Experimental RSFs of 163 Dy (left panel) and 164 Dy (right 
panel). The dashed lines in the uppermost panels show the contri- 
butions from the giant dipole resonances multiplied by K, the solid 
line represents the sum of these two contributions, k(Je\ +/mi)- 
The fit to the experimental datapoints in 164 Dy is performed up to 



(11) = 5.3 MeV. In the lower panels the sum k(/ei + /mi) is sub- 



tracted from the experimental data and the fit to the Ml pygmy reso- 
nance is displayed (solid line). 



is displayed. We notice that the fit to the experimental data 
around the Ml pygmy resonance is good, especially for 
163 Dy. When compa ring the pygmy resonance parameters of 
i63,i64 Dy (see Table |in| to those extrac ted for i60,i6U62 Dy 

reported in [7], we find a smaller width of the pygmy 
resonance. The previous measurements for 16 °. 161162 Dy 



yielded widths in the range of Tpy = 1 .26 — 1 .57 MeV using 
a constant temperature of Tf = 0.3 MeV. In the present 
work, with the same constant temperature, we find widths 
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FIG. 5: [Color online] The integrated B(Ml f) strength of the pygmy 
resonance for dysprosium isotopes measured with the Oslo method 
(filled circles) and their average value (solid line). The TSC data 
point for 163 Dy (8| is also displayed (filled square). 




0.80 MeV for 163 Dy and 



of Tpy = 0.86 MeV and T py 
164 Dy, respectively. The nucleus 163 Dy has been investigated 
earlier by the Prague group, analyzing TSC spectra from the 
162 Dy(ra,2y) 163 Dy reaction 0. In their work, the width of 

: 0.6 MeV. For 



the pygmy resonance was reported to be r py 
this specific case ( 163 Dy), the measured F 



Y energy E (MeV) 



FIG. 6: Experimental RSF for 160 Dy. The dashed line represents the 
tail of the GEDR, while the dashed-dotted lines give the contributions 
from the GMDR and the Ml pygmy resonance. The solid line is the 
sum of all the resonances. The fit to the experimental datapoints is 



py 



from the Oslo performed up to E y = 6.9 MeV. 



data and the data from the Prague group is comparable within 
the uncertainties. 

We note from Fig. [4] that Omi for 163 Dy is significantly 
larger than for 164 Dy. The reason for this is not yet under- 
stood. In order to obtain a more precise comparison, the total 
integrated strength B(M1 f) given by 



B(Ml t) 



9hc 



32n 2 \ E J 



(12) 



Ml pygmy 



is calculated for 16 °- 164 TJy, and the results are displayed in 
Fig. [5] When calculating the weighted average of the Oslo 
data, a value of 6.6(4) ^ is found 1 . The B(M1 f) value from 
the TSC experiment is not included in the fit because no errors 
are given in Ref. [8j. We observe that all the measured values 
agree within the uncertainties. 

For 164 Dy, we observe an increase in the RSF compared 
to theory for energies above Ey ss 5.0 MeV. Similar features 
have been observed in (7, /) experiments on other nuclei, e.g., 
116124 Sn (ED and 208 Pb (23 . For these nuclei the structure is 
thought to be governed by the so-called neutron skin oscil- 
lation, a collective mode of El character that for stable nu- 
clei is located in the region of Ey = 5 — 10 MeV. This feature 
has been observed in nuclei with a high neutron-to-proton ra- 
tio N/Z, and is interpreted as an oscillation of the neutron- 
enriched periphery of the nucleus versus a core consisting of 



equally many protons and neutrons, N = Z |2][26l. Enhanced 
strength is also observed in the RSF of 117 Sn measured at 
OCL ll27ll . Unfortunately, the present experimental technique 
cannot provide information on the electromagnetic character 
of the enhanced strength in 164 Dy. However, it might be a 
reasonable guess that the observed strength stems from the 
E 1 skin oscillation, since we note that Dy nuclei have a high 
neutron-to-proton ratio of N/Z = 1.36 — 1.48 for the stable 
isotopes. Evidence of both the Ml pygmy resonance and the 
E 1 pygmy resonance in one and the same nucleus has, how- 
ever, not been reported earlier. 

Data on 160 Dy from a previous experiment [7] also appears 
to have excess strength, see Fig. [6] Unfortunately, the strength 
function in the interesting region (Ey > 7 MeV) suffers from 
poor statistics. However, this could be a hint that the same 
feature is present in this nucleus. 



IV. SUMMARY AND CONCLUSIONS 

The nuclei 163 > 164 Dy have been investigated using the Oslo 
method. The radiative strength functions have been extracted, 
displaying the Ml pygmy resonance. This resonance has been 
studied in detail, and it is found that the measured widths are 
smaller than what has previously been measured in other Dy 



nuclei at OCL. However, the pygmy widths of 



163,164 



Dy in 



1 The B(M1 f) values for l«>.l6l,l62 Dy are ca l cu l at ed from Ref. (6). The 
values of 161162 Dy are the weighted averages of the values obtained from 



the ( 3 He,a) and ( 3 He, 3 He') reactions. 



the present work are still larger than what has been measured 
for 163 Dy by means of the TSC method. When comparing the 
total integrated strength B(M1 f) of the Ml pygmy resonance, 
the results for all the nuclei agree within the uncertainties. 
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For 164 Dy, we have observed an excess of strength for 
Ey > 5 MeV compared to model calculations; similar fea- 
tures can also be seen in 160 Dy. The enhanced strength might 
be due to neutron skin oscillations. If that is the case, this 
is the first time both the scissors mode and the neutron-skin 
oscillation mode is seen in one and the same nucleus. 
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